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Simulating Simplified Versions of the IKKT Matrix Model 
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We simulate a supersymmetric matrix model obtained from dimensional reduction of 4d SU(A'') super Yang- 
Mills theory (a 4d counter part of the IKKT model or IIB matrix model). The eigenvalue distribution determines 
the space structure. The measurement of Wilson loop correlators reveals a universal large N scaling. Eguchi- 
Kawai equivalence may hold in a finite range of scale, which is also true for the bosonic case. We finally report 
on simulations of a low energy approximation of the lOd IKKT model, where we omit the phase of the Pfaffian 
and look for evidence for a spontaneous Lorentz symmetry breaking. 



1. Motivation 

Several candidates for a constructive definition 
of superstring theory have recently attracted at- 
tention. Here we focus on the IKKT model (or 
IIB matrix model) [|l|, which is supposed to cor- 
respond to IIB superstring theory in the large N 
limit. That correspondence is supported by some 
analytical arguments. This matrix model is for- 
mally obtained from ordinary super YM gauge 
theory in the zero volume limit (one point). 

Here we study directly the large N dynamics 
of large N reduced matrix models. Some re- 
sults were obtained before for the "bosonic case" 
(where the fermions are dropped by hand), but 
we now want to address mainly the SUSY case. 
In particular we simulate the 4d counterpart of 
the IKKT model |§ — we denote it as 4d IIB 
matrix model — which corresponds again to the 
dimensional reduction of 4d super YM gauge the- 
ory. This model has also been studied analytically 
|3| and numerically in the framework of dynami- 
cal triangularization [g] . Here we report on direct 
Monte Carlo simulations using the Hybrid R al- 
gorithm 0. Conclusive results can be obtained 
because in the 4d version the fermion determinant 
is real positive — in contrast to the lOd IKKT 
model, where simulations would be plagued by a 
sign problem. We got away with a computational 
effort of 0{N^) in the SUSY case, and of 0{N^) 



in the bosonic case. 

2. The 4d IIB matrix model 

The 4d IIB matrix model is given by 



Z = 

Sb = 
Sf = 



dA e'^" / d^d4> e-^' 

-^Tr(v;ar^^K,#]) 



(1) 



where A^, ip, 
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V'a (a^ = 1 ■ • ■ 4, a = 1, 2) are com- 
plex, traceless N x N matrices, and the A^j, (only) 
are Hermitean. We use T — ia, r4 — H. In ad- 
dition to SUSY and SO (4) invariance, this model 
has a SU(A^) symmetry, which is inherited from 
gauge invariance. 

The first question about this model is if it is 
well-defined as it stands. Since the integration 
domain of dA is non-compact, divergences are 
conceivable. However, our results, as well as re- 
sults on a number of special cases [^-g] confirm 
consistently that this model is well-defined at any 
N; there is no need to impose a IR cutoff. This 
implies that the only parameter g is simply a 
scale parameter, rather than a coupling constant. 
It can be absorbed by introducing dimensionless 
quantities. The challenge is, however, to tune g as 
a function of N so that the correlators are finite 
at N ^ oo, see Sections 4 and 5. 



3. The space structure 

In the IIB matrix model, the space coordinates 
arise dynamicaUy from the eigenvalues of the ma- 
trices Afj^ |l] . In general the latter cannot be di- 
agonalized simultaneously, which implies that we 
deal with a non-classical space. We measure its 
uncertainty by 



N 



Tr(A2) 



7nax 
UgSU{N) 



]{iuA,u^ur 



and the "maximizing" matrix U is also used for 
introducing the coordinates of N points, 

iUma.A^Ul,,)u il = l...N). (2) 



^l,fl 



What we are really interested in is their pairwise 
separation r{xi, Xj) — \xi — Xj\, and we show the 
distribution p{r) in Fig. 1. We observe p ~ at 
short distances {r/y/g< 1.5), hence a UV cutoff is 
generated dynamically. We also see that increas- 
ing N favors larger values of r. To quantify this 
effect we measure the "extent of space" [| 



^new 



rp{r) dr 



(3) 



Fig. 2 shows Rnew and A as functions of N (at 
g — 1). The inclusion of fermions enhances Rnew 
and suppresses A, keeping their product approxi- 
mately constant. The lines show that both quan- 
tities follow the same power law, Rnew, A ex 
iVi/-*, in SUSY and in the bosonic case. In SUSY 
this behavior is consistent with the branched 
polymer picture: there one would relate the num- 
ber of points as A^ ^ {Rnew/(-Y" , where (. is some 
minimal bond, which corresponds to the above 
UV cutoff. The Hausdorff dimension d/f = 4 then 
reveals consistency with our result. 

4. Polyakov and Wilson loops 

We define the Polyakov loop P and the Wilson 
loop W — which is conjectured to correspond to 
the string creation operator — as 



P{p) 
W{p) 



— Trfe'P^' 

TV 



(4) 



iTrf 

TV V 



gipAigipA2g-ipAig-ipA2 \ t^\ 



•^Note that the quantity B? = I r"^ p{r) dr- diverges, 
since p{r) oc r~^ at large r (see second Ref. in |6|). 



P(r) 




N=48 


0.014 


s^^^ 


32 


0.012 


-it 'M 


24 


0.01 
0.008 
0.006 


16 




o^/ %^% 




0.004 


1* % 




0.002 


- f ^^ 







J . , ^ 





5 10 15 20 r/rg 

Figure 1. The distribution of distances between 
space-points in the SUSY case at various N . 
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Figure 2. The "extent of space" Rnew cind the 
space uncertainty A as functions of N at g = 1 . 
Of course the choice of the components of v4^ is 
irrelevant, and the parameter p G IR can be con- 
sidered as a "momentum" . 

Now g{N) has to be tuned so that (P), (W) 
remain finite as iV — > cx). This is achieved by 

(6) 



g ex l/VN 



which leads to a beautiful large N scaling; Fig. 
3 shows the invariance of (P) for A^ = 16 ... 48 
in SUSY. Also the bosonic case scales accurately 
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Figure 3. The Polyakov function in the SUSY 
case for various values of N and g-'N = const. 

The historic 2d Eguchi-Kawai model had a 



{Zn)^ symmetry, which impUed {P{p ^ 0)) = 0, 
a property which was crucial for the proof of the 
Eguchi-Kawai equivalence to gauge theory. As we 
see, this property is not fulfilled here, but {P{p)) 
falls off rapidly, towards a regime where the as- 
sumption of this proof holds approximately. 

We proceed to a more explicit test of Eguchi- 
Kawai equivalence by checking the area law for 
{W{p)). Fig. 4 shows that the area law seems 
to hold in a finite range of scale. Remarkably, 
the behavior is very similar B] in the bosonic 
case. There we further investigated the behavior 
at much larger N |l^, and we observed that the 
power law regime does neither shrink to zero — as 
it was generally expected — nor extend to infin- 
ity — a scenario which seems possible from Fig. 
4. At least in the bosonic case its range remains 
finite at large N. Recently the same behavior was 
observed in a study of the lOd bosonic case |11|. 
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Figure 4. The Wilson loop in the SUSY case for 
various values of N and g^N — const. 

5. Multipoint functions 

We now consider connected multipoint func- 
tions (©102 ■ ■ ■On)con, Oi being a Polyakov or a 
Wilson loop. We wonder if it is possible to renor- 
malize all of those multipoint functions simply by 
ZOi, so that a single factor Z 



[ren) 



(ren) ^r\{ren) 



.o. 



(ren)\ 



inserting Ol 

renders all functions {Of^^'O. 

(simultaneously) finite at large N. 

It turns out that such a universal renormaliza- 
tion factor seems to exist in SUSY. We have to 
set again g c>c 1/vN, and then Z oc N provides 
large N scaling, as we observed for a set of 2, 3 
and 4-point functions. Two examples are shown 
in Fig. 5. Our observation can be summarized by 
the SUSY rule 

(0)=0(1), {Oi...O,,) = 0{N-") (n>2). 



This implies that large N factorization holds, 
(Oi...O„) = (Oi)...(0„> + 0(^"'), as in 
gauge theory, although coupling expansions are 
not applicable here. 




Figure 5. A Polyakov 2-point function 

(C'f^ = ([ImP(fc)]2); and a Wilson 3-point 
function (Cf^ = {[linW {k)fY{.eW [k)) - 
{\hnW{k)f) (ReVK(fc))J (both connected), with 
g^N = const, and renormalization factor Z oa N , 
which leads to large N scaling in the SUSY case. 

For the bosonic case, a 1/d expansion H sug- 
gests large N factorization to hold as well, but it 
also predicts (d ...On) = 0(7V-2("-i)) (^ > 
2). This is confirmed numerically 0: in partic- 
ular the 3-point functions shown in Fig. 5 now 
requires Z^ ex N'^. Therefore no universal renor- 
malization factor Z exists in the bosonic case, 
which is an important qualitative difference from 
the SUSY case. 

6. Simulations in ten dimensions 

We also performed simulations in d = 10 |12[ , 
where we simplified the IKKT model as follows: 

(1) We use a 1-loop approximation, which is ex- 
pected to capture the low energy dynamics. This 
amounts to an effective action, keeping track of 
off-diagonal elements only to the quadratic order, 
in the spirit of Ref . ||] . 

(2) We omit the phase of the Pfaffian by hand, 
in order to avoid the sign problem. Thus a Monte 
Carlo study becomes feasible. 

The validity of (1) is supported by our results 
for Rnew, but (2) is certainly a drastic step. Still 
one could hope to observe basic properties of the 
IKKT model at least qualitatively. These sim- 
plifications allow for a simulation effort of only 



Our main interest here is if the eigenvalue dis- 
tribution of the A^ indicates a spontaneous sym- 
metry breaking (SSB) of SO(IO) invariance, as 
it was suggested in the formulation of the IKKT 
model 0. 

To this end, we consider the moment of inertia 



T„, 



N{N^ 



-T \ / A'^'M Xjfi)[Xw Xji/) (I) 



i>j 



{i = 1 . . . A^), and we measure its 10 eigenvalues. 
A gap in this spectrum would indicate the SSB of 
Lorentz symmetry. However, this cannot be ob- 
served, even though we raised N up to 512. On 
the contrary, we observe that the eigenvalue dis- 
tribution becomes more and more isotropic as N 
increases, see Fig. 6, and the same is true in d = 6 
|l^ . We conclude that if SSB of Lorentz symme- 
try occurs in the IKKT model, then it must be 
driven by the imaginary part of the action. 
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Figure 6. The spectrum of the moment of inertia 
(normalized by VN) for N = 192, 256, 384, 512. 
7. Conclusions 

We first simulated the 4d IIB matrix model, 
both, SUSY and bosonic. We varied N up to 
48, which turned out to be sufficient to study the 
large N dynamics. 

We confirmed that the model is well-defined 
as it stands, hence g is a pure scale parameter. 
The space coordinates arise from eigenvalues of 
the bosonic matrices A^t. The extent of space 
follows a power law. In SUSY this agrees with 
the branched polymer picture. Fermions leave 
the power unchanged but reduce the space un- 
certainty — though it remains finite at large N. 

The large N scaling of Polyakov and Wilson 
loops and their correlators requires g oc 1/\N 
in SUSY and in the bosonic case, but the wave 



function renormalization is qualitatively different: 
only in SUSY a universal renormalization exists. 

The area law for Wilson loops holds in a finite 
range of scale for the SUSY and the bosonic case. 
The latter comes as a surprise, and we checked up 
to rather large N that this range remains indeed 
finite. Hence Eguchi-Kawai equivalence to ordi- 
nary gauge theory H may hold in some regime. 

Finally we simulated a lOd low energy effective 
theory, where the phase was dropped by hand. 
We could not observe any sign of a spontaneous 
breaking of Lorentz symmetry. 

W.B. would like to thank for the excellent organiza- 
tion of LATTICE 2000, where he enjoyed discussions 
with A. Gonzalez-Arroyo, U.-J. Wiese and J. Wosiek. 
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